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Abstract 

The phonon thermal contribution to the melting temperature of nano-particles is inspected. 
Unlike in periodic boundary condition, under a general boundary condition the integration volume 
of low energy phonon for a nano-particle is more complex. We estimate the size-dependent melting 
temperature through the phase shift of the low energy phonon mode acquired by its scattering on 
boundary surface. A nano-particle can have either a rising or a decreasing melting temperature 
due to the boundary condition effect, and we found that an upper melting temperature bound 
exists for a nano-particle in various environments. Moreover, the melting temperature under a 
fixed boundary condition sets this upper bound. 
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I. INTRODUCTION 


This work focuses on the quantum size effect with particular emphasis on environment- 
depending melting temperature of nano-particles. The continuing progress in the design of 
nano-particles led to enhanced and novel functionality^. Concerning the thermal stability of 
heat resistant 2 , we ask a question: can nano-particles stand a temperature higher than the 
one its fixed surface partners melt at? There are many different thermodynamic theories of 
small systems^ - -, each considering different important aspects of size-dependent melting of 
nano-particles. The discrete quantum energy level has not been carefully considered in these 
theories. The critical role of phonon in thermal related phenomena is well-known^’ 7 . But 
only recently the finite spacing energy levels became a greater awareness in size-dependent 
melting^. There is a chance to understand the rising of size-dependent melting temperature 
of nano-particles with more attention on the important aspect of quantum finite size effect. 
We estimate in this work the change of melting temperature by coating a fixed size particle 
or changing its environment. We give the upper limit of the change on melting temperature. 

The melting temperature for small particles was modeled and studied in more than 100 
years ago^, and the pressing need for a deep understanding continues today^~— . In this 
study, we use the Lindemann melting criterion^ - — for its simplicity in estimating the melting 
temperature of a nano-particle with different boundary conditions. By Lindemann criterion, 
a nano-particle melts at the temperature T m at which the ratio of u, the square root of the 
mean square of atom thermal displacement, to a the lattice constant reaches the Lindemann 
critical value L c : 

_ u(T m ) 

L'c 

a 

II. PERIODIC BOUNDARY CONDITION 

We repeat some necessary derivation for u{T) in harmonic approximation. A lattice 
specified by a set of the vectors Ri will be studied, with one atom at each lattice point. The 
displacement u, of an atom from its equilibrium position R, can be calculated by 

H = — U ia + 2 Rj )Ui a UjQi (2) 
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where a = x,y,z and Ui a is the ath component of the displacement. Ri = (i x - 2 ^v ~ 

-gd, i z — ^rp)a, i Q = 1 The center of the nano-particle is at the origin of the 

coordinates, and the boundary of the nano-particle is set on R a = ±4 surfaces in coordinate 
space. N = L 3 is the number of atoms of the nano-particle. M is the mass of the atom. The 
potential energy <f> is expanded to the second order and < h Q / 3 (Ri, Rj) = (d 2 §/du ia dujp) 0 . 

The equations of motion of the lattice is then 


Mui a ^ ^ ^^(Rj, Rji^yn 
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The general solution can be written in vibration modes Qka- 
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ca CT (k) and eta are phonon frequency and phonon polarization of wave-vector k, respectively. 
The atomic mean-squared thermal displacement is^ 
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where () means grand canonical ensemble average. The center-of-mass motion is u cm = 
The square root of the mean square of atom thermal displacement, u{T ), is given 


by : 


u(T) = 
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Melting temperature T m is obtained by solving Eq. mi. We will use boundary conditions to 
account for various environments. 

For periodic boundary condition the mode expansion function for atom displacement in 
each a-direction is f(k a ,R ia ) = s J^e lkaRia , k Q = 2n Q 7r/La, and n a = —L/ 2 + 1, • • •, L/2. 
The melting temperature T mn for size L nano-particle under a periodic boundary condition 
is then calculated^: 
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where u;(k) = ay^k) is used when summing up the three polarized vibration directions in 
each k mode. u cm is given only by the zero wave-vector phonon at k — 0. The second 
term removes contribution of the global moving ( u cm ) 2 . This is the missing of phonon 
contribution from the zero-mode volume. We have repeated briefly the derivation of this 
zero-mode volume in the work of Sui et al^. For the bulk system, melting temperature T mb 
is the solution of the same equation in the limit of 1/L — > 0, which turns the second term 
into zero: 

2 _ 3 ha r /a 3 1 

27t 3 M J 0 ^(kjtanht^ggy 

III. BOUNDARY REFLECTIVE PHASE SHIFT 

We analyze the phase shift of the low energy phonon mode for different boundary con¬ 
ditions. A theoretical estimated expression for the melting temperature of nano-particle 
will contain physical quantities with some uncertainty. Beside the main variable, the size 
L, other physical factors may play a role in finite size melting, such as shape, surface re¬ 
construction, and environmental effect etc. An expression has to have physical quantities 
not so certain to absorb all the different factors from material to material, environment to 
environment, and particle to particle. We will use the phase shift of low energy acoustic 
phonon for this function. A model with boundary is good for nano-particles under different 
boundary conditions. The phase shift of phonon mode depends on its scattering on the 
surface of a nano-particle. A phase shift can be drawn from the experimental data, or it has 
to be chosen with some uncertainty for a particular situation before we can make estimation 
for a design of the melting temperature. 

The scattering and phase shift have been studied most thoroughly in the context of 
quantum theory. The notion of reflection of waves in one-dimensional scattering plays a 
central role in the following detailed discussion. We model the boundary of the nano-particle 
as an additional potential which is bigger than zero outside the nano-particle and zero inside 
the particle. In any a-direction, the general solution for Eq. Q is in terms of waves moving 
in opposite directions, f(k a , Rio) = c e lkoiRia + d e~ lkaRia , for each vibration mode Qka- We 
study the standing waves under a general real boundary potential. The most general form 
for the reflection coefficient would then be r = e zSka ; the termination at the boundary could 
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at most introduce a phase change in the reflected wave: 


^ g ikafta. 
Q gikotRoi 



R a = La/ 2. 


( 9 ) 


Different boundary conditions resulting in the same magnitude of phase shift are physi¬ 
cally equivalent. A very careful but long and technical analysis on phase shift and boundary 
condition can provide considerable understanding of nano-particle’s surface and environ¬ 
ment. We will put this study under control. We assume an isotropic surface and environ¬ 
ment. The same phase shift will be assumed in each a-direction. The phase shift for a small 
wave-vector is expanded up to the first order of 1/L: 5k a = 5 + a\k a . 5 is a constant phase 
shift and cq is the expansion coefficient for the first order term. We will discuss on phase 
shift in the range of [—7r, 0], with corresponding boundary barrier effectively repulsive. This 
phase shift is the parameter we used to model the boundary effect. 

The wave-vector is fixed by the boundary condition Eq. (|9]). First, when the standing 
wave is of even parity, we have c/d = 1 in f{k a) Ri a ), and e 1< ' kaLa+Ska ' ) = 1. The mode 
expansion for atom displacement on this kind of phonon mode is f(k^R ia ) ~ cos (k^R ia ), 
with k e a = 2n ^+\^a\ ' qq ie q^eger n a runs from 0 to L/2 — 1. The second group is the 
odd parity group for c/d = —1 in f(k a ,Ri a ), and e^ kocLa+n+Ska ^ = 1. The mode expansion 
for atom displacement on this kind of phonon mode is /(/c°,i?,: Q ) ~ sinwith k° a = 
(2n a -i)TT+\5 k o a \ ' integer n a runs from 1 to L/2. The low energy wave-vector has an 

increase \Sk a \/La when <5 fcct ^ 0. 


IV. FIXED BOUNDARY CONDITION 

With the above inspection, we will be able to estimate the melting temperature of a nano¬ 
particle under a general boundary condition. The low energy wave-vector has an increase 
\8ko\/La. The melting temperature is increased by a boundary condition shifting k a wave- 
vectors upward. The discrete summation of phonon modes in Eq. (J5]) can be written into 
I( 2 \s\-Tr)/ 2 La^ka for each cc-component of wave-vector-. When |5fc Q | = tt/2 the integration 
volume for low energy phonon starts at k a = 0: f Q dk a , which results in the same amount 
of low energy phonon contribution to atom displacement as the case for bulk material in 
Eq. (JHJ). The cases for |5 fca | < tt/2 will depress the melting temperature of a nano-particle. 
The lower bound of this depression when |5*. a | —> 0 was carefully studied in the work of Sui 
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et. al^. When |<5fc a | > 7 t/ 2, the melting temperature of a nano-particle will be higher than 
its bulk parent. 

We study the rising of the melting temperature due to the phase shift. The boundary 
scattering shifts wave-vectors up to k a = (|5 fea | + n a 7r)/La by a nonzero Sy. a . The range of 
the phase shift Sk a in k a is [— 7T, 0]. If |(5*, a | = 7t/2, the density of the low energy phonon 
states in a nano-particle is the same as the density in the bulk material. We prefer a smaller 
thermal displacement of the atom when we study a heat resistant. We want less phonon 
contribution to atom displacement. The least contribution from phonon modes occurs when 
bka = —K- This situation is given by the ideal fixed boundary condition (fixed-BC), for 
which the boundary sets up a node surface for standing acoustic waves. 

The upper melting temperature bound is set by the ideal fixed-BC for nano-particles. The 
phonon wave-vectors move up for this boundary condition and some more low energy phonon 
contribution is taken away by the deduction of the second ( u cm ) 2 term. In the picture of 
microscopic physics, less low energy phonon will contribute to atom displacements in finite 
size fixed-BC particle than in bulk material. Therefore, the melting temperature increases for 
the boundary effect and the finite size effect. Moreover, for all possible boundary conditions 
modeled in phase shifts, the low energy phonon modes shift up the most for fixed-BC, and 
their contribution missed the most for fixed-BC. Thus the melting temperature increases 
the most for fixed-BC. 

The upper limit for effective phase shift is |<5fc Q | = w, and the melting temperature in¬ 
creases at most when a particle subject to a fixed boundary condition. However, this ideal 
boundary condition is not easy to achieve when preparing nano-powder heat resistant, since 
there is no infinitely heavy walls we can coat microscopically on a nano-particle. Fixed-BC 
is an unrealistic boundary condition but it is instructive. Nano-particles can be grown inside 
sieves, can be coated with other elements. The melting temperature of a nano-particle can 
approach but not reach the one for fixed-BC. From the above discussion on both finite size 
and boundary condition, we can conclude that the upper melting temperature bound is set 
by the ideal fixed-BC for nano-particles. 

We derive this upper melting temperature bound in the following. For fixed-BC, the 
mode expansion for atom displacement is f(k^,R ia ) = cos(k^R ia ) for even parity, and 
f{k° a ,R ia ) = \J^si\\{k 0 a Ri a ) for odd parity. Quite different from periodic condition’s mode 
space, there is no k — 0 phonon mode for fixed-BC. Even parity will contribute to 
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u cm , since | Y^ia=i cos[k^a(i a — ^p]| = 1/sin The melting temperature T mn for size L 
nano-particle under a fixed-BC is then calculated from Eq. (JT]) : 


3 h 



( 10 ) 
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( k ) tanhf^^] Ha l Ls in ^] 2 


The first term is the summation of uf a defined by Eq. ([5]), in which k a = n a 7i / La with 
n a = 1, • • ■, L. The second term removes the contribution of the global moving (u cm ) 2 , 
which sums only the even parity k e a = (2 n a + l)7r /La with n a = 0, • • •, L/ 2 — 1. 

V. DISCUSSION AND APPLICATION 

State-of-art functional designs of nano-particles developed in recent years^ 1 ^. The first 
general factor came into our understanding on size-dependent melting was the surface volume 
ratio of number of atoms? 1 . Nano-materials are studied case by case in experiments and Ab 
initio calculations^ - —. With the help of computer simulation, in principle all the parameters 
can be fixed to certain extent for a particular design of a particular particle, of its surface 
and its properties'!^. A designed melting temperature could be accurate up to 1 K—. 
Mechanical stable region for nano-materials is also under current studying along with the 
estimation of T mn , and giving the safe working temperature^. Many on going theoretical 
researches are working for more general understanding on size-dependent melting. 

Statistical mechanics is valid for an ensemble of big systems and small systems. Thermo¬ 
dynamics of small systems is well understood. New and powerful expressions were derived 
through thermodynamics^ for size- and boundary-dependent melting, but they have not 
converged to physically equivalent ones yet. Size-effect on intensive variables can be in¬ 
cluded through thermodynamic equations or ensemble average of quantum states. In this 
study on general boundary conditions for nano-particles, it is clear that the boundary con¬ 
dition and size together determine the wave-vectors of phonon modes, and so determine 
the low energy density of states of phonon. Such a size- and boundary-dependent density 
of states was not included in previous approximations in a phonon Debye temperature. In 
real applications, a particle interacts with other surrounding atoms, molecules, and particles 
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Fig. 1 fixed-BC 



FIG. 1: The size dependent melting temperature for nano-particles under fixed-BC. The lines and 
points are T mn /T m b with and without the second term in Eq. (fTOl) . Guiding lines simply connecting 
the points are presented. 

simultaneously. Our analysis is right when the bulk phonon dispersion is almost kept in a 
nano-particle^. The individuality of the particle is retained under the statistical average. 
Vibration modes inside the particle see only an averaged reflective phase shift. 

In this study, a simple procedure to estimate the melting temperature for nano-particles 
under a general boundary condition can be suggested. In this simple estimation procedure 
we use sound speed v to approximate the phonon dispersion for all phonon modes in Eq. (jHJ) 
for an integration, and in Eq. (TTOjl for a discrete summation. For example, in Fig. 1 we 
plot the size-dependent melting for a simple cubic lattice under fixed-BC. Phonon dispersion 
is approximated by o;(k) = sin 2 + sin 2 ^ + sin 2 We substitute hv/ak B T mb by 
0.25 when solving Eq. (jSJ) and Eq. (fTOl) for T mn /T mb at each size L. 

Fig. 1 shows that the global movement contributes more in the variation of melting 
temperature as the size decreases down to L = 20 and smaller. The global movement 
deduction accounts for five percent of the change when L = 10. For a nano-particle at the 
size of L = 20, with hvjoksT-mb ~ 0.25 for its parent bulk material, we can read on Fig. 
1 that an increase of melting temperature less than 18 percent can be expected when we 
design a heat resistant coating for the nano-particle. This example shows that the reflective 
phase shift at the particle boundary plays an important role in size-dependent melting, and 
size and phase shift are both general factors in this phenomenon. 
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VI. CONCLUSIONS 


We investigated in this paper the melting temperatures for nano-particles under a general 
boundary condition. We find that the missed contribution from the low energy phonon will 
raise the melting temperature. This phenomenon comes from finite size effect, and boundary 
condition plays an important role in it. The melting temperature for the fixed boundary 
condition gives the upper bound for the melting temperature of a nano-particle in all kinds 
of environment. This upper bound is not reachable in any real nano-particle design. 
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